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Abstract 

A calculation of the spin-wave polarization operator is very important for the 
analysis of the magnetic structure of high temperature superconductors. We 
analyze the significance of the incoherent part of the spin-wave polarization 
operator within the framework of the t-J model. This part is calculated 
analytically for small doping with logarithmic accuracy. We conclude that 
the incoherent part of the spin-wave polarization operator is negligible in 
comparison with the coherent part. 
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I. INTRODUCTION 

It is widely believed that the spin structure of the ground state of high-temperature cop- 
per oxide superconductors is of the type of the Anderson spin-liquid state On the other 
hand it is well known that at zero doping the materials under consideration are insulators 
with long-range antiferromagnetic (AF) order. The instability of long-range AF order under 
doping in the t-J model was probably first realized by Shraiman and Siggia 0. It was 
later demonstrated in numerous works (see e.g. refs. The instability is due to the 

strong interaction of spin-waves with mobile holes. The spin-wave polarization operator is 
needed for the analysis of this instability and especially for the derivation of the spin-liquid 
ground state. The calculation of the so called coherent part of the polarization operator is 
straightforward. However there is also the incoherent part. The purpose of the present work 
is to calculate this part and to provide a rigorous proof of its smallness in comparison with 
the coherent part. It means that in the analysis of the magnetic properties one can neglect 
the incoherent contribution to the spin-wave polarization operator. This conclusion is in 
agreement with that of ref. H] based on the numerical computation. 



II. HAMILTONIAN 

The t-J model is defined by the following Hamiltonian: 

H = H, + Hj = -t Yl {dLdma + ii-C.) + J E Sn-S„, (1) 

<nm>a <nm> 

where dl^^ is the creation operator of a hole with spin a (cr =t, i) at site n on a two- 
dimensional square lattice. The d]^^ operators act in the Hilbert space with no double 



electron occupancy. The spin operator is S„ = \dl^^(Tai3dnp. < nm > are pairs of nearest- 
neighbor sites on the lattice. Below we set J as well as the lattice spacing equal to unity. 

At half-filling (one hole per site) the t-J model is equivalent to the Heisenberg AF model 
P,pi]|] which has long-range AF order in the ground state |TT,T^. Let us denote the wave 



function of this ground state by |0). This is the undoped system. We consider the doped 
system based on the ground state of the undoped system. In spite of the destruction of the 
long-range AF order it is convenient to use |0) and corresponding quasiparticle excitations 
as a basis set for the doped system. The effective Hamiltonian for the t-J model was derived 



in terms of these quasiparticles in the papers [13-15 



^eff = ek/^L^ka + ^q("i<^q + /^J/^q) + Hh,sw + Hhh- (2) 



ko- 



It is expressed in terms of normal spin-waves on an AF background ttq, /3q (see e.g. review 
|T^) and composite hole operators /iko- (o" = =tl/2)- The summations over k and q are 
restricted inside the Brillouin zone of one sublattice where 7q = ^{cosqx + cos qy) > 0. The 
spin-wave dispersion relation is 



^q = 2^1 - 72 , uj^ ^ v^|q| for q < 1. (3) 



The properties of single holes are well established (for a review see ref. IjT^])- The wave 
function of a single hole can be represented as ipica = ^ko-|0)- large t the composite hole 
operator hl_^ has a complex structure. For example at t/J = 3 the weight of the bare hole 
in T/^kcr is about 25%, the weight of configurations of the type "bare hole -|- 1 magnon" is 
~50% and for configurations of the type "bare hole -|- 2 or more magnons" it is ~25%. 
The dressed hole is a normal fermion. The hole energy has minima at k = kg, where 
ko = (±7r/2, ib7r/2). For t < 5 the dispersion can be well approximated by the expression 



18 



ek ^ ^0 + 2 - + 4.56t2 - 2.8t272 + -/52(cos K - cos kyf. (4) 

The numerical values in this formula are some combinations of the Heisenberg model spin 
correlators. The constant Eq defines a reference level for the energy. To find Ck with respect 
to the undoped system one has to set -Eq = 0. However for the present work it is convenient 
to set eko = and therefore Eq = ^0.66^ -|- 4.56^^ — 2. The coefficient (32 is small and 
therefore the dispersion is almost degenerate along the face of the magnetic Brillouin zone 



7k = 0. According to refs. [|19| , |20| ~ 0.1 1 for t > 0.33. To avoid misunderstanding we 
note that formula (^) is not valid for very large t where the hole band width is saturated 
at the value of the order of unity and does not increase with t. However, physically we are 
interested in t «i 3 (see e.g. refs. pT|-[23[1) where (HI) works well. Near the band minima, kg 



the dispersion (H) can be presented in the usual quadratic form: 

1.0 1 



ep ~ ^(3ipi + P2<^(3u (5) 

where pi (^2) is the projection of k — ko on the direction orthogonal (parallel) to the face 
of the magnetic Brillouin zone (fig. 1). From eq. (^ for t 3> 0.33 we find (3i ~ 0.65t, hence 
the mass anisotropy is (3i/(32 ~ 7. For small hole concentrations, i.e. 5 -C 1, the holes are 
localized in momentum space in the vicinity of the minima of the band ko = (±7r/2, ±7r/2) 
and the Fermi surface consists of ellipses. (We recall that 6 = corresponds to an insulator 
with long-range AF order.) The Fermi energy and Fermi momentum of non- interacting holes 
are 



and Pf ~ y/PiFP2F ~ (vr(5) 



1/2 



(6) 



The Fermi momentum pp is measured from the center of the corresponding elhpse. 

The interaction of a composite hole with spin- waves is of the form (see e.g. refs. P^J2^JTB[ ]^ 



Hh,sw = Zl5'k,q(4+qi^kTaq + 4+qT^ki/5q + H.C. 



k,q 



(7) 



5fk,q = 2V2/(7kt/q + 7k+qVq) 2^^^ f ^{q^ sin + Qy sin ky) for g < 1. 



Here [/„ 



7^ + i and Va 



sign(7q),/^j — i are the parameters of the Bogoliubov 



transformation diagonalizing the spin-wave Hamiltonian. The hole spin-wave coupling con- 



stant / is a function of t evaluated in the work [|13[. For large t the coupling constant is 
t-independent, / ~ 2. Let us stress that even for t > J the quasihole-spin-wave interaction 
(|^ has the same form as for t <^ J (i.e. as for bare hole operators) with an added renormal- 
ization factor (of the order J/t for t ^ J). This remarkable property of the t-J model is due 
to the absence of a single loop correction to the hole-spin-wave vertex. It was first stated 



implicitly by Kane, Lee and Read In refs. [jT9|,^,|T3[ it was explicitly demonstrated that 
the vertex corrections with different kinematic structure are of the order of a few percent 
at t/J ^ 3. There is also some g-dependence of the coupling constant /. For example 
fiq 



L15/(g = 0) at t/J = 3 (see refs. |[T5|j2q| ). However this dependence is weak 



and is beyond the accuracy of the calculation of the renormalized value of /. Therefore we 
neglect this dependence. 

Finally there is the contact hole- hole interaction Hhh in the effective Hamiltonian (Q). Hhh 



is discussed in detail in refs. |ll^J15| , |26[| . It is proportional to some function A(t). For small 
t this function approaches the value —0.25, which gives the well known hole-hole attraction 
induced by the reduction of the number of missing antiferromagnetic links. However for 
realistic superconductors t ~ 3 (see e.g. refs. pl|-p3|). Surprisingly the function A(t) 
vanishes exactly at t ~ 3 and this means that the mechanism of contact hole-hole attraction 
is switched off. In contrast the spin-wave exchange mechanism H^ g^ is negligible for small 
t for which f ^ t and it is most important for large t for which / approaches 2. We are 
interested in "physical" values of t: t ~ 3. Therefore in the present work we neglect the 
contact interaction {H^h = 0) and consider only the hole-spin-wave interaction at t = 3. 
The corresponding value of / according to ||13[ is / = 1.8. Thus all numerical values in the 
present work are calculated at t = 3, /3i//32 = 7 and / = 1.8. 



III. POLARIZATION OPERATOR 



There are two types of spin-waves in the t-J model, a and f3. Therefore, one generally has 
to introduce a set of spin-wave Green's functions |^^: Daa, Dap, D^a, and -D/3/3. However 
if we restrict ourselves to the long wavelength limit, the consideration can be simplified. It 
is well known that in the long wavelength limit the Heisenberg model is equivalent to the 
nonlinear a- model (see e.g. review paper |]l6l). Therefore the usual field theory crossing 
symmetry is valid. The same is valid for the t-J model. This is evident from eq. (|^): at 
g -C 1 the vertex g^ ^ ^ 5'k-q,q- This means that instead of considering a set of Green's 



functions one can introduce a single combined Green's function of vector excitation: 



where Il{uj,q) (or P{uj,q)) is the mobile holes polarization operator. For the system to be 
stable the condition 



cjq + 2n(0, q) > (9) 
should be fulfilled. Otherwise the Green's function (|) would possess poles with imaginary 

(V. 

In the single loop approximation for dressed quasiholes the polarization operator is given 
by the diagram in fig. 2. The quasihole Green's function in this diagram is determined by 
the effective Hamiltonian (0). Let us call this contribution to the polarization operator the 
"coherent part" . The meaning of this term is elucidated below. In this paper we neglect 
the interaction between quasiholes and consider them as a normal Fermi-liquid. In this 
approximation the calculation of a single loop polarization operator is very simple. In the 
static long wavelength limit the result is [0 

neo^(0,q)^ ^-y^q for g « p^. (10) 

This is independent of the hole concentration! Using this result one can easily prove that 
the stability condition (^ is violated because —211(0, q)/co'q = 2/^/(7rA/3^^) ^ 2.8. This is 
a well known fact concerning the instability of long-range AF order under doping (see e.g. 
refs. M). 

Let us now explain the obtained result in terms of bare dna which enter into the bare 
Hamiltonian (|l|). The Green's function of a bare hole is of the form 

G(e, k) = + Gincoh. (11) 

e - Ck + zO 

where Z is the quasiparticle residue. If we use the bare hole Green's function we also 
have to use the bare hole-spin- wave coupling constant /tare = 2t in the vertex (7k,q (|3). 
Substitution of the pole (coherent) part of (|TT]) into the single loop polarization operator 
gives the combination fbareZ, but this is exactly the effective coupling constant / [|r3|. Thus 



in the single loop approximation the effective theory with Hamiltonian (0) only takes into 
account the coherent (pole) part of the bare hole Green's function (|TI|) . Because of this the 



contribution ( [T0| ) is called the coherent part. 

The incoherent part of the polarization operator which comes from Gmcoh in (|Tl|) is 
proportional to the hole concentration 6. Therefore it is negligible for 5 ^ L This conclusion 
agrees with that of Becker and Muschelknautz . 

Now we will calculate Pmcoh explicitly. We give two different calculations of Pincoh- The 
first calculation is based on the self-energy approach. It is well known that the self-energy as 
well as the polarization operator is in essence an energy shift of the system. Let us consider 
the diagram in fig. 3a which represents the energy shift of a hole (solid line) in an external 
spin- wave (wavy line). If we join up the wavy lines and integrate over the momenta of the 
spin-wave we get the usual hole self energy (fig. 3b) which is already taken into account in 
the correct (renormalized) hole dispersion (H) and the renormalized value of the spin-wave 
coupling constant /. If we join up the solid lines and integrate over the momenta of the 
hole we get the coherent polarization operator (fig. 3c) which has also already been taken 
into account. Let us consider now the first correction to the diagram fig. 3a. It is given by 
the diagrams in fig. 4. (We remind the reader that there is no one loop correction to the 
spin- wave vertex.) This is the correction to the hole energy in an external spin- wave or the 
correction to the energy of the spin-wave in the presence of a single hole. The calculation 



of this correction is very simple if one uses Schrodinger perturbation theory. Consider for 
example the case when the external spin-wave is of the /3-type, i.e. for this wave the z- 
projection of spin is 5^ = +1. The arrows in fig. 4 indicate the projections of hole spin (note 
that there is a spin flip at each vertex). It is evident that for spin up diagram 4b vanishes 
and similarly for spin down diagram 4a vanishes. We will see that the integral over Q is 
logarithmically divergent. Therefore we set g <^ Q <^ 1. Near the band minimum the hole 
energy is quadratic in momenta (^. Therefore it can be neglected in comparison with the 
energy of the intermediate spin-wave cjq which is linear in momentum. This means that the 
energies of all three intermediate states in fig. 4a as well as in fig. 4b are approximately 
equal to cuq. Therefore the energy correction corresponding to fig. 4 is independent of the 
hole spin projection and equals 

The integration over Q is restricted by the hmits Qmin ^ Pf = Vtt5 and Qmax ~ tt. 
The integral over Q is logarithmically divergent and therefore it gives a big logarithm: 
\ia{Q max /Qmin) ~ 0.51n(l/5). When the momentum of the hole is close to the band mini- 
mum, k ko = (7r/2, ±7r/2), one gets 

Se=J^^-^^^lnS. (13) 
4\/27r q 

To find the incoherent spin-wave polarization operator we now need to sum (p!3| ) over all 
occupied hole states (i.e. sum over k). Thus for u eF and q <^ Pf the result is 

P^ncoh{^,c^)^ -^^■6-\n--q. (14) 

There are also higher order diagrams which contribute to Pincoh- An example is presented 
in fig. 5. However all of these diagrams do not contain the big logarithm (ln(l/5)) and 
therefore can be neglected according to the accepted accuracy. 

To elucidate the meaning of Pincoh we will now give a different derivation of (|I^) . The 
first correction to the single loop polarization operator arises from the diagram presented in 
fig. 6. After integration over the energies in the closed loops one gets the following expression 
for the contribution of this diagram: 

/ \ f ^ ^ ^ 2 2 / \ 

^^^^'^^ = 7 {2n)H2ny '^'^'^'-'' ^^^^ 

^k-q(l - ^k-q) - ^k(l - ^k-q) ^ ^k^q^k-q - ^k^k-q ^ 

ek-q - Ck + u^)^(ek-q - Ck-q -ujq + uj) (ek-q - ek + u;)^(ek-^q - Ck-q + ujq + uj) 

-nk(l - T^k-q) ^ (1 - nk)rak-q 

(ck-q - ek + tuq)^(ek-q - Ck-q + ujq - uj) (ck-q - Ck - t^q)^(ek-q - Ck-q - uq - u 

where rik is the ground state occupation number. The first two terms in the expression 
(|T5|) have second powers of (ck-q — e]^ + u) in the denominator. This factor is exactly the 
denominator of the one-loop polarization operator fig. 2. It means that the first two terms 
correspond to the renormalization of the pole (coherent) part of the hole Green's function 
(pJ]). They are actually already taken into account in Pcoh because for the calculation of 
Pcoh we have used the dressed hole Green's function which absorbs all of the corrections to 



the pole part. Thus the new contribution, corresponding to Pincoh, arises from the last two 
terms in (|r^). After carrying out some algebraic manipulation and using the fact that q is 
small (and so, for example ek-q ~ ek) Pincoh can be written as 

P..co.(0,q) = -2j (2vr)2(27r)2^'^'Q^'^''»(u;Q + ek-Q-6k)3- ^^^^ 

It is evident that this expression is equivalent to (|12|). Neglecting Ck-q — Ck in comparison 
with cuq and carrying out the integration over Q with logarithmic accuracy we once more 
come to the result (111) for Pincoh- 



IV. CONCLUSION 

In the present work we have analytically calculated the incoherent contribution to the 
spin-wave polarization operator. This analytical calculation is possible because of the loga- 
rithmic enhancement of the corresponding diagrams. The ratio of the incoherent contribution 
(|l^) to the coherent one (0) is 



mcoh(0, q) 1 „2 



P m ^ ^/V/^i/?2-5-ln- . (17) 
Pcoh{^,fV 8 " d 

For the parameters of the t-J model corresponding to realistic high-T^ superconductors and 
for hole concentration 5 < 0.2 this ratio does not exceed 10%. It means that in the analysis 
of the magnetic properties one can neglect the incoherent contribution to the spin-wave 
polarization operator. 
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FIGURE CAPTIONS 

FIG. 1. The magnetic Brillouin zone. 

FIG. 2. The coherent spin-wave polarization operator. 

Fig. 3. a) The leading order contribution to the energy shift of a hole (sohd hne) in 
an external spin-wave (dashed hne). b) The hole self-energy. c)The coherent spin- wave 
polarization operator. 

FIG. 4. The fourth order correction to the energy of a hole in an external spin-wave. 
FIG. 5. A typical fifth order diagram for the energy of a hole in an external spin-wave. 
FIG. 6. The two loop correction to the spin-wave polarization operator. 
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